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We consider an effective Lagrangian containing contributions from glueball and gluon degrees of 
freedom with a scale-invariant coupling between the two. The thermodynamic potential is calculated 
taking into account thermal fluctuations of both fields. The glueball mean field dominates at low 
temperature, while the high temperature phase is governed by low-mass gluon-like excitations. The 
model shows some similarities to the lattice results in the pure glue sector of QCD. In particular, 
it exhibits a strong first order phase transition at a critical temperature of approximately 265 MeV 
' when reasonable parameters are taken. 

O ' 
(N 

£ 

i— » 

CO 



- 1—1 

X 



I. INTRODUCTION 



The aim of this paper is to non-perturbatively describe the pure glue sector of quantum chromodynamics (QCD) with 
an effective Lagrangian containing a small number of parameters. In this endeavor we are guided by the requirement 
I ' that the model reproduce the most important non-perturbative features of QCD. At zero temperature the model 
should therefore account for the gluon condensate. As the temperature increases through a critical temperature, T Cl 
{SJ • a phase transition should occur in which the condensate is melted and gluon-like excitations become the relevant 
t-H , degrees of freedom. This picture is supported by lattice calculations which are quite well established for the pure 
OO ■ glue sector Above T c the thermodynamic "data" has been fitted by several simple models ranging from one 

OS ■ employing massless gluons with a low-momentum cut-off || to a model assuming temperature dependence in both 
the effective gluon mass and the bag constant [Q (see this reference for a more complete review). 

In this paper we formulate a simple effective field-theoretical model for constituent gluons which become massive 
via an interaction with the gluon condensate. The constituent gluons are described by an Abelian vector field, while 
the gluon condensate is identified with the dynamical glueball field which has a non-zero expectation value. The form 
of the glueball potential is well known Q , fixed by the requirement that scale invariance be broken through generation 
of a trace anomaly as in QCD. Therefore it is surmised that the non-linearities of the actual QCD Lagrangian can 
be modeled by the glueball potential and the mass-like coupling term. Below we demonstrate that this simple model 
can indeed reproduce some of the features of QCD. 

The plan of this paper is as follows. In Sec. II the effective Lagrangian is written down and the fields are 
decomposed into mean field and thermally fluctuating parts. Expressions are given for the equations of motion, 
masses and thermodynamic variables. Since we have to deal with a non-linear potential we introduce a novel method 
of handling the fluctuations. In Sec. Ill we discuss the choice of model parameters and present our results. Sec. IV 
is reserved for conclusions and outlook. 



II. THE MODEL 



A. Effective Lagrangian 

Motivated by the above considerations we write our effective Lagrangian as 

C = \d^cj> - U(<j>) - ±A M „ • A"" + iG 2 2 A M • A^ , (1) 

where c/> is the scalar glueball field and A M is the Abelian vector field with strength tensor A M „ = d^A. u — d„A^. The 
dimension of the vector A^ is a priori unknown, and we define it to be v/3 so that the effective number of constituent 
gluon degrees of freedom is v. The second term in Eq. (fil) is the glueball potential 
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where A is a constant and A, which is the only dimcnsionful parameter in the model, defines the vacuum glueball field, 
0o = A/e^. The last term in ([l]) gives a masslike coupling between the <f> and A M fields with a coupling constant G. 

Since G is dimensionless, scale invariance is broken only by the glueball potential, U(4>). This potential was 
constructed to reproduce the trace anomaly of QCD in an effective theory Indeed, for our Lagrangian the trace 
of the energy-momentum tensor is 

<% = 4U-4^ = -\*4>\ (3) 

so that the vacuum energy density is — \\ 2 <j>o = — 2 -Bo- In QCD, the latter is proportional to the gluon condensate and 
is roughly known from QCD sum rules ||; see Narison Q for a more recent value. Therefore a nonzero expectation 
value for the glueball field (cf>) implies the presence of a gluon condensate. In the following we consider ((f)) to be the 
order parameter for our study of the thermodynamics and phase structure at finite temperature. The model thus 
contains only four parameters: Bo, </>o and G. The pure glueball sector of the Lagrangian ([!]), which we will treat 
as a special case, has been discussed previously by Agasyan J9|. 

The Lagrangian ([!]) is of similar structure to the dual Ginzburg-Landau model |8| . The latter also assumes Abelian 
dominance, for which lattice calculations offer some support, and considers three scalar magnetic monopole fields, 
rather than the single glueball field here. Where a quartic potential has been assumed in Ref. g], we take a logarithmic 
potential for reasons previously explained. In both cases the square of the constituent gluon mass is proportional to 
2 , or equivalently the square root of the gluon condensate, a dependence common in other models as well 



B. Equations of Motion and Masses 

It is convenient to take the ratio of the glueball field to its vacuum value, \ — "fif^o- Then the Euler-Lagrange 
equations of motion are 

d 2 A ll - (d v A v ) + g 2 X 2 K = , (4) 

where we have defined g = G(f>o. The non-trivial point of our treatment is that in addition to the mean field we include 
the thermal fluctuations of the glueball and gluon fields in a consistent way. To that end we break the glueball into 
mean field and fluctuating parts, x — X + A with (A) = 0, where the angle brackets denote a thermal average. In the 
second of Eqs. (^) we replace x 2 by its thermal average (x 2 ) so that we can interpret the third term as a mass term; 
this amounts to imposing the condition d" 'A„ = 0, as appropriate for a vector field with three degrees of freedom. 

The field fluctuations are decomposed into plane waves and we take the thermal average of the equations of motion, 
assuming that the thermal average of the product of glueball and gluon fields can be approximated by the product of 
their respective thermal averages. This gives the mean-field versions of Eqs. (^) for the glueball field, 

2£o(x 3 logX 2 ) = 5 2 x(A M -A^), (5) 

and the subsequently vanishing mean gluon field, (A M ) = 0. 

The dispersion relations for the gluon and glueball excitations, including contributions from thermal fluctuations, 
can be written 

e\ = k 2 + m*l ■ e 2 x = k 2 + m* 2 , (6) 
where the effective masses are defined to be the thermal average of the second derivative of the potential, i.e. 

<t>tm* x 2 - - = SB {x 2 Inx 2 ) + 4B (x 2 + (A 2 )) - g 2 (A, ■ A") 

< = (|^)-5 2 (X 2 + (A 2 )). (7) 

This set of equations is closed by expressing the quantities (A^-A^) and (A 2 ) in terms of the field quanta distributions. 
Using standard methods one obtains for the vector and scalar fields, respectively, 
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[A ll -A>') = ~fdk^n B {e A ) ; (A 2 ) = ^ J dk ^«M< x ) • (M 
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where tib(x) — [eP x — is the Bose-Einstein distribution function and (3 = 1/T is the inverse temperature. The 
equation for the mean field (|B|) and the equations for the masses ([?]) must be solved self-consistently. 



C. Evaluation of the Thermal Fluctuations 

Equations (|^) and (^) require the thermal average of functions of x which involve a logarithm and these are 
handled in the following manner. Consider a general function f(x) and Taylor expand the fluctuations so (f(x)} — 
So° / < ""' ) (x)(A ,1 )/?t.!, where /W denotes the n th derivative of the function. Here we are considering the contributions 
from a single vertex, ignoring loop diagrams with two or more vertices, as is appropriate for a mean field treatment. 
Simply truncating the Taylor expansion at low order as in Ref. S would be appropriate for low temperatures, but 
would be inadequate for high temperatures where the mean field vanishes and the fluctuations are large. Fortunately 
we can treat the problem exactly. Taking the thermal average of each possible pair of fields A for a given rt-point 
vertex gives (A") = (n — 1)!!(A 2 )^ for n even and zero for n odd (see Ref. jll| for further discussion). Defining a 
Gaussian weighting function 

P(z) = (2^<A 2 )pexp-(^y , (9) 

one finds that 

(f(x)}= I <!- /'C- I I > /""U'-T 1 : / + HOI 
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where in the last step we have resummed the Taylor series. Thus the fluctuations enter with a Gaussian weighting 
and Eq. ( |l0| ) is straightforward to compute for any \. By performing the series expansion one can see that x = is 
an exact solution of Eq. (||). In this case the integral (|lo|) can be carried out analytically |Q and Eqs. (0) become 

2 m* 2 =6i?o(A 2 )ln a <A 2 }- 5 2 <A M .A^ ; m^ 2 = 5 2 <A 2 } , (11) 
where In a = I — 7 — In 2, with Euler's constant denoted by 7, giving a ~ 4.0402. 



D. Thermodynamics 

The grand canonical potential per unit volume can be written in a straightforward way: 

00 

£ = iSo(x 4 (lnx 2 - !)> + \B - ^ 2 <(A 2 ) - JL| rffcfc^^A) + ns^y {u) 



Here a constant second term has been added so that fl = at zero temperature and the third term subtracted so as 
to avoid double counting. Notice that, apart from the constant term, all the quantities are temperature dependent, 
for instance (x 4 (mx 2 — |)) is evaluated using Eq. (|Io| ) which involves the temperature-dependent quantity (A 2 ) of 
Eq. (|). 

In order to have consistent thermodynamics f2 must be a minimum with respect to variations in the mean field \- 
Performing the minimization we indeed find that Eq. (|^) is the necessary condition. Thus our equation of motion, 
our mass equations and the grand potential treat the thermal fluctuations in a coherent and well defined manner. 
The pressure is simply P = —Cl/V, and the energy density is easily obtained: 

00 

+ 27r 2 J dkk2 ( veArlB ( eA ) + e x n B(e x )) ■ ( 13 ) 


Note that in the case \ = 0, (x 4 lnx 2 ) = 3(A 2 ) 2 In a(A 2 ), where a is the constant defined previously. 
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III. RESULTS 



A. Choice of Parameters 



The model has four free parameters: v, Bq, cf>o, and G. The effective number of gluon degrees of freedom, v, 
determines the asymptotic behaviour of the equation of state. In order to have £/T 4 rs 4.7 at high temperature, as 
found on the lattice [|, we need to choose v w 14 for gluons. The standard degeneracy for massless gluons in SU(3) 
is, of course, 16. We will also consider v = 6 which would roughly correspond to SU(2), as well as the case where 
gluons are excluded and we have a pure glueball theory [y — 0). For a given v the critical temperature, T c , of the 
phase transition (see below) is largely determined by the quantity -Bo- For SU(3) we take v = 14 and choose a value 
for Bo of (391 MeV) 4 so as to reproduce the deconfinement temperature found in lattice calculations. However Bq 
also determines the zero-temperature vacuum gluon condensate and our value must be consistent with independent 
determinations of this quantity. Our Bq is somewhat larger than the old value of (340 MeV) 4 found by Shifman, 
Vainshtein and Zakharov ||, but in good agreement with the more recent updated average of (399 ± 13 MeV) 4 given 
by Narison Q] . The magnitude of the vacuum energy density associated with our value of the gluon condensate (bag 
constant) jBq w 0.8 GeVfm -3 . The third parameter, </>q, can be fixed by appealing to the vacuum glueball mass, 
m x = 4i?o/0o> which follows from Eq.(7|) at T = 0. For the glueball mass, values in the range 1.5 - 1.7 GeV are 
suggested by data and by calculations |l3]; for definiteness we take m x =1.7 GeV following Sexton et al. [Q. Finally 
in order to fix G we assume that the glueball is a loosely bound system of two gluons and therefore choose the effective 
gluon mass in vacuum, m,A = G(f>o = <?, to be \m x \ the recent study of phenomenological gluon propagators |l5| ] 
suggests that this is a reasonable estimate. 



B. Phase Transition 



Fig. 1 shows the mean glueball field, x, for the three cases mentioned above. Since x an d the other variables are 
essentially constant at lower temperatures the abscissa starts at T = 100 MeV. A most striking feature of the model 
is that it exhibits a first order phase transition at a critical temperature T c . Here the mean glueball field drops from 
a value of slightly less than unity to zero, and it remains zero for T > T c . The arrows on the figure indicate where 
the transition takes place (the remainder of the curves correspond to metastable or unstable branches). The physics 
behind this is indicated in Fig. 2, where we plot il/V (relative to the value at \ = 0) as a function of x for various 
temperatures in the vicinity of T c . Below T c there exists a local minimum of the effective potential at \ — 0, but the 
absolute minimum is at x ~ 1- As the temperature is increased the latter minimum becomes shallower and at T c it 
has the same depth as the minimum at % = 0. Thereafter the stable solution corresponds to x = and ultimately 
the minimum at x ~ 1 disappears. 

For the pure glueball case (dotted line in Fig. 1) the transition temperature, T c = 490 MeV, is much higher than 
that estimated by Agasyan ||] . This is probably due to the fact that in that work the potential was expanded to low 
orders, whereas here we have an essentially exact treatment of the thermal fluctuations. When the glue degrees of 
freedom are included the transition temperature is much reduced and, using reasonable values for the parameters, is 
found to lie in the neighborhood of the SU(3) lattice result T c — 270 ± 5 MeV for the pure glue sector of QCD [ {To) . 
For our chosen parameters the first order transition for the solid curve occurs at 265 MeV. For v = 6 the dashed curve 
shows that the critical temperature rises to 340 MeV. To compare with SU(2) one should also take into account the 
scaling of the gluon condensate with the number of colors. Reducing the condensate by a factor of | results in a value 
of T c = 300 MeV. These figures are reasonable in view of lattice calculations which yield a 20% increase in T c (with 
considerable error) |l7| when the number of colors is changed from three to two. However the lattice results indicate 
a second order transition for SU(2) [Q, whereas we have a first order transition, although mean field treatments are 
expected to be inadequate in the neighborhood of critical points. 



C. Effective Masses 



The ratios of the effective masses to the temperature, m x /T and m* A /T, are displayed in Fig. 3 (the complicated 
structure of the low temperature metastable region for v — 14 is not relevant here and for clarity is suppressed in the 
figures). The masses change little below the critical temperature and so the ratio drops as the temperature increases. 
Beyond the critical point the masses grow linearly with temperature to a good approximation. Such behavior is 
expected at very high temperatures where perturbation theory is applicable. For the solid curves (y — 14) m* A /T 
is quite small, 0.18, whereas m*/T is large, 4.5, so that the glueball plays only a minor role beyond the transition 
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temperature, as is physically expected. The reason that one inevitably has a large glueball mass follows from Eq. 
(fTT|). The gluon contribution to to* (second term) is positive and large at high degeneracy v, so in order to have 
a low glueball mass the first term would have to be negative. This would be only be possible if the argument of 
the logarithm were less than 1. However for a low-mass glueball with to*/T <C 1 the argument of the logarithm 
a(A 2 ) ~ [T/(1.72(/>o)] 2 . So in order to have a negative value for the logarithm for T > T c , (j>o would have to be 
substantially increased compared to the chosen value of <po — 180 MeV. This would lead to a reduction of the vacuum 
glueball mass to an unreasonably low value. 



D. Thermodynamics 

The pressure and the energy density are shown in Fig. 4, where we plot £/T A and 3P/T 4 . At low temperatures 
the energy density and pressure are little changed from the vacuum values. The phase transition temperature is 
determined by the point at which the pressure curves for the two stable solutions intersect. We note that the model 
predicts interesting behaviour for the pressure in the vicinity of this critical point. On cooling from high temperatures 
the system can enter a supercooled metastable phase, which can even have zero pressure. It is conceivable that in 
relativistic heavy ion collisions such behavior could allow metastable, supercooled droplets of hot gluonic fluid to be 
produced fig|j . 

Beyond T c the pressure and energy density in Fig. 4 increase rapidly to their asymptotic values. This qualitative 
behavior for the pressure is in agreement with lattice calculations, although the approach to the asymptotic value is 
sharper in our model. At high temperatures 2>P « £ and the value of these quantities for the solid and dashed curves is 
very close to that expected for an ideal gas of massless gluons with the appropriate degeneracy, v — 14 or 6, indicating 
that we have reached the asymptotic regime. One aspect of this model is seemingly in strong disagreement with the 
lattice data. Namely, at the critical temperature there is a large latent heat and the energy density overshoots the 
asymptotic value of £/T 4 , approaching it from above. This is in contrast to the lattice calculations where the latent 
heat is a factor of 2-3 smaller |jl],|o) and £/T i approaches its asymptotic value from below. 

For T > T c the behavior can be understood rather simply. The ratio of the glueball mass to temperature, to*/T, is 
large, while the gluon mass is small and, to a good approximation can be taken to be zero. This allows the necessary 
Bose integrals to be approximated easily (see, for example, Ref. f2l}| ). It turns out that only the constant term and 
the thermal gluon term are numerically significant for the pressure and energy density. Thus, one obtains expressions 
very similar to those of the bag model 

3P _ tt 2 v 3B £ ^ tt 2 i/ B 

T 4 ~ ~30 4T 4 ' T 4 ~^0 _ + 4T 4 ' ^ ' 

Since the pressure is approximately zero at the phase transition, the above equation provides an estimate of the 
critical temperature T c « {45So/(27r 2 i/)} 1 / 4 which is accurate to 10% or better. The thermodynamics of the model 
given by Eq. (|lj) can be contrasted with a power law fit to the lattice results Jl| at high temperatures not too close 
to T c , which gives 

3P\ . „ „ /T c \ 216 / £ \ , „ _ /T, N 



Indeed the sign of the deviation of the energy density from the asymptotic value is opposite to the prediction of 
the bag model. This results in too large a latent heat, as we have mentioned earlier. We remark that the dual 
Ginsburg-Landau model appears to suffer from similar difficulties. 

The qualitative features of our model are not sensitive to variation of the parameters within reasonable limits. For 
example, the primary effect of varying the glueball-gluon coupling constant G is to scale the gluon effective mass 
according to Eq. (Q). There is little change in T c , and above this temperature the thermodynamics is still represented 
by Eq. ( J14| ) since m* A /T remains small (unless G is made unreasonably large). Similarly, the main effect of altering <f>o 
is to change the masses; increasing <j>a decreases the glueball mass and increases the gluon mass. This assumes that the 
constant B has been fixed to vacuum expectations. If we choose to alter B , it is the vacuum (and low-temperature) 
glueball mass which must change. This alters the critical temperature according to T c ~ B^ 4 , as we have mentioned. 
Above T c the value of Bq is inconsequential for the masses, while the thermodynamics still follows Eq. (|l4|). 

We have also considered introducing an (A M • A p ) 2 term in the Lagrangian, which would be suggested by QCD, but 
this does not appear to improve the situation. However it is worth noting that the simple replacement to*^ 2 — > to^+to 2 ,, 
where too is a constant mass of order 500 MeV, can alter the predictions. With some modest adjustment of the other 
parameters, the energy density and pressure above T c can be brought into semi-quantitative agreement with the 
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lattice data. In particular, the energy density no longer overshoots the asymptotic value. This modification, however, 
is inconsistent with QCD in that the corresponding term in the Lagrangian, ^TOqA^ • A M , introduces a dimensional 
parameter and consequentially gives an unwanted contribution to the trace anomaly. Thus such an addition contradicts 
our initial approach. It does however strongly suggest that the problem with our model is that above T c the masses 
immediately become proportional to the temperature and some important physics has not been accounted for in our 
description of the masses in the region T c < T < 2T C . 



IV. CONCLUSIONS 



In conclusion, we have examined a very simple Lagrangian model for the pure glue sector of QCD and determined 
its thermal properties using a novel treatment of the thermal fluctuations. In accord with physical expectations, the 
model shows a phase transition between a low-temperature, glueball-dominated regime and a high-temperature phase 
dominated by low-mass gluon-like excitations. Sensible parameters yield a transition temperature in agreement with 
lattice simulations. Beyond T c the qualitative behavior of the pressure is reasonable, but the energy density in the 
neighborhood of the phase transition disagrees with the lattice data. This pathology appears to be a result of the 
the low effective gluon mass immediately after the phase transition. So while the thermodynamics are reasonable at 
low and high temperatures, the lattice results reveal our model to be incomplete just above the critical temperature. 
Nevertheless, we conclude that our simple Lagrangian indeed captures important non-perturbative features of QCD. 

In the future it would be interesting to consider modifications suggested by the color dielectric model which simulates 
confinement. It would also be worthwhile to consider the inclusion of quark degrees of freedom. Such effective models 
may be particularly useful for simulations of time-dependent processes and the effects of finite baryon density which 
are yet intractable on the lattice. 
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FIG. 1. The mean glueball field, x, as a function of temperature, T. For the dotted curve the gluon field is neglected, while 
the dashed and solid curves correspond to degeneracies, v, of 6 and 14 for the gluon field, respectively. The arrows indicate 
where a phase transition takes place and the thermodynamically stable phase becomes x — 0. 
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FIG. 2. The grand potential density, [SI — Sl(x = 0)]/V, as a function of the mean glueball field, x, for temperatures, T, in 
the vicinity of the critical temperature T c . 
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FIG. 3. The effective masses of the glueball and gluon in units of the temperature, m*/T and m* A /T, as a function of 
temperature. See caption to Fig. 1 for the meaning of the curves. 
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FIG. 4. The thermodynamic quantities £/T 4 and 3P/T 4 as a function of temperature. See caption to Fig. 1 for the 
meaning of the curves. 
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